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Abstract 


Two tropical squall lines from TOGA COARE and GATE were simulated 
using a two-dimensional cloud-resolving model to examine the impact of 
surface fluxes on tropical squall line development and associated precipitation 
processes. The important question of how CAPE in clear and cloudy areas is 
maintained in the tropics is also investigated. Although the cloud structure and 
precipitation intensity are different between the TOGA COARE and GATE squall 
line cases, the effects of the surface fluxes on the amount of rainfall and on the 
cloud development processes are quite similar. The simulated total surface 
rainfall amount in the runs without surface fluxes is about 67% of the rainfall 
simulated with surface fluxes. 

The area where surface fluxes originated was categorized into clear and 
cloudy regions according to whether there was cloud in the vertical column. The 
model results indicated that the surface fluxes from the large clear air 
environment are the dominant moisture source for tropical squall line 
development even though the surface fluxes in the cloud region display a large 
peak. The high-energy air from the boundary layer in the clear area is what feeds 
the convection while the CAPE is removed by the convection. The surface 
rainfall was only reduced 8 to 9% percent in the simulations without surface 
fluxes in the cloudy region. Trajectory and water budget analysis also indicated 
that most moisture (92%) was from the boundary layer of the clear air 
environment. 
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1 . Introduction 


It is well known that sensible heat, latent heat, and momentum fluxes 
between the ocean and the atmosphere play an important role in cloud 
development and precipitation processes over the ocean. Surface fluxes are 
temporally and spatially complex in the region of active convection. For 
example, earlier observational studies (Malkus and Ronne, 1954; Malkus, 1958; 
LeMone and Pennell 1976) have demonstrated that formation of oceanic 
cumulus clouds in the trade winds is controlled largely by marine boundary layer 
processes. Observational studies in the Western pacific warm pool region 
(Bradley, 1991; Young et al., 1992; Fairall et al., 1996) have shown that surface heat 
and momentum fluxes all have a peak in the convective leading edge due to 
strong gust winds and colder air temperatures in the convective region. The 
surface fluxes in the large clear area are much smaller and more uniform than 
those in the convective region. Several numerical modeling studies (Tao et al., 
1991; Wang et al., 1996) have indicated that sensible and latent heat fluxes can 
enhance surface precipitation and cloud coverage by comparing simulations 
with and without the effects of ocean fluxes for both subtropical (TAMEX) and 
tropical (TOGA COARE) squall lines. Wang et al. (1996) also showed that among 
the heat and momentum fluxes, the latent heat flux is the most important 
component for cloud development. However, those studies left following 
question unanswered, which part of the surface fluxes, those in the large clear air 
environment or those in the convective area, is responsible for enhancing cloud 
development and precipitation? 

In this sensitivity study, simulations were made for two well documented 
tropical squall lines, the 12 September 1974 GATE (Szoke and Zipser, 1986) and 
the 22 February 1993 TOGA COARE cases (Jorgensen et al., 1997). Both cases, 
TOGA COARE and GATE, (Table 1) have moderate convective available 
potential energy (CAPE), 1400 and 1600 J/kg, respectively. Tropical oceanic 
convective systems are typically associated with a moderate CAPE. While the 
TOGA COARE case has a very moist environment with a precipitable water of 
6.05 g cm“2 , the GATE case is substantially drier with a precipitable water of 4.80 g 
cm‘2. The sea surface temperature in the TOGA COARE case is higher that that 
in the GATE case. The surface heat fluxes are strongly dependent on the 
temperature difference between the air and the sea sureface. The environmental 
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winds are also quite different between the two cases. In the TOGA COARE case, a 
fairly strong low level jet (about 12 m s'l) is present at a height of 2 km, and there 
is a weak overturning upper tropospheric wind (4 ms ') at about 10 km. The 
GATE case has less shear in the lower troposphere, but there is a strong jet in the 
upper troposphere above 10 km (about 30 m s _1 ) in the same direction as the low 
level flow. 

The objective of this study is to identify which part of the surface fluxes, the 
clear large-scale environment or the cloudy region with strong gust winds, is 
more influential in the development and precipitation processes of tropical 
squall lines. The important question of how CAPE in clear and cloudy areas is 
maintained in the tropics is also investigated. The two well-documented tropical 
squall lines cases from GATE and TOGA COARE, that occurred in quite different 
large scale environments, are simulated and the results are analyzed. 

2. Model and Model Setup 

The tool used in this study is the two-dimensional version of the Goddard 
Cumulus Ensemble (GCE) model. The simulated flow is assumed to be 
anelastically balanced. Sound waves have been filtered out by neglecting the local 
variation of air density with time in the mass equation. The cloud microphysics 
include a parameterized Kessler-type two-category liquid water scheme (cloud 
water and rain), and a parameterized Lin et al. (1983) or Rutledge and Hobbs 
(1984) three-category ice-phase scheme (cloud ice, snow and hail/ graupel). Short- 
wave (solar) and long-wave (infrared) radiation parameterizations (Chou 1984; 
1986) as well as a subgrid-scale turbulence (one-and-a-half order) scheme are also 
included in the model. 

The GCE model has implemented a Multi-dimensional Positive Definite 
Advection Transport Algorithm (Smolarkiewicz and Grabowski, 1990). All 
scalar variables (potential temperature, water vapor, turbulence coefficient and 
all five hydrometeor classes) used forward time differencing and the MPDATA 
for advection. The dynamic variables, u, v and w, used a second-order accurate 
advection scheme and a leapfrog time integration (kinetic energy semi- 
conserving method). A stretched vertical coordinate (height increments from 40 
to 1150 m) with 31 grid points was used to maximize resolution in the lowest 
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levels of the model. A total of 1024 grid points were used in the horizontal with 
750 m resolution. Details of the model description can be found in Tao et al. 
(1993), Tao and Simpson (1993), and Simpson and Tao (1993). 

The surface flux parameterization used in this study is from the TOGA 
COARE flux algorithm (Fairall et al., 1996). This is primarily based on the bulk 
scheme developed by Liu et al. (1979), which has shown good agreement with 
observations (Bradley et al, 1991). The transfer coefficients for momentum, 
sensible heat, and latent heat fluxes are based on the Monin-Obukhov similarity 
theory of the atmospheric surface layer (Businger et al, 1971). This bulk scheme 
has been modified (Fairall et al, 1996) to better simulate surface fluxes from the 
tropical ocean and calibrated by the large amount of data from the TOGA 
COARE field experiment. 

The shortwave radiation models of Chou (1992) are used to compute solar 
heating in the atmosphere and clouds and at the surface. The longwave radiation 
of Chou and Suarez (1994) is used to compute cloud and atmospheric infrared 
cooling. The cloud optical properties are parameterized using a broadband 
emissivity method (Stephens, 1984). Both the liquid and solid phase of the water 
as well as their size distribution are used in the parameterization. Details of the 
cloud optical calculations and their sensitivity tests, as well as a review of cloud 
resolving modeling studies on the cloud-radiation interactions can be found in 
Tao et al (1996). 

The soundings used for this sensitivity study are from two well- 
documented squall line cases, the 22 February 1993 TOGA COARE (Wang et al., 
1996) case, and the 12 September 1974 GATE case (Ferrier et al., 1995). The set-up 
for the model domains and boundary conditions for both cases are the same. 
The model domain consists of 1024 horizontal and 31 vertical grid points with 
open lateral boundary conditions. A constant horizontal grid spacing of 0.75 km 
in the interior 864 grid points is placed within a coarser, horizontally stretched 
outer region with a stretching ratio of 1.05, resulting in 18506 km of total 
horizontal model domain with 684 km being constant resolution. The vertical 
coordinate is stretched in order to maximize resolution in the lowest levels of 
the domain. The vertical grid spacing increases from 40 m near the surface to 
1164 m at the top of the 21.5-km deep domain. 
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3. Sensitivity test of surface fluxes 


The experiments designed to test the influence of surface fluxes on tropical 
squall lines from different regions is listed in Table 2 (T for TOGA COARE and G 
for GATE). The control runs (T1 and Gl) have all the model physics, including 
surface fluxes over the whole computational domain. Runs T2 and G2 have no 
surface fluxes anywhere in the entire simulation domain. The surface fluxes 
were set to zero in the cloudy region in runs T3 and G3, and the surface fluxes i n 
the clear region were set to zero in the runs T4 and G4. The clear and cloudy 
regions were determined at every time step using a total hydrometeor content of 
10' 5 g/g as a threshold value. Every column in the domain was searched 
vertically to determine whether the column was clear or cloudy. If it exceeded 
the threshold at any level, it was considered cloudy. The largest amount of cloud 
coverage at any time in any run was about 1.6% (300 km) of the simulation 
domain including both TOGA COARE and GATE cases. 

3.1 Surface precipitation characteristics 

Squall line convection can alter the sea surface fluxes in the cloudy area. Fig. 
1 shows the instantaneous surface momentum and heat fluxes at 10 hours into 
the GATE control simulation (Gl). Similar flux characteristics for the TOGA 
COARE case were shown in Wang et al., (1996). The average latent heat flux in 
the clear area for TOGA COARE (80 wm’ 2 ) is greater than that in the GATE case 
(60 wm-2). This is primarily due to the fact that the TOGA COARE case has a 
larger air-sea surface temperature difference. Both cases show a large jump in 
fluxes at the leading edge of the squall line systems due to the strong surface 
winds. Following the peak, there is a transitional area where the surface fluxes 
gradually decrease to the values of the environment. 

The domain averaged total surface rainfall amounts from the various 
simulations are listed in Table 3. The TOGA COARE and GATE runs give 
consistent results: the least rainfall was from the runs without surface fluxes, and 
the most rainfall was from the control cases. Although the magnitude of the 
fluxes in the convective region is large due to stronger surface winds, the surface 
fluxes from the clear area have a much greater influence on the surface rainfall 
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than the fluxes from the cloudy area. The runs without surface fluxes in the 
cloudy region have a 5 to 10% rainfall reduction, verses a 26% rainfall reduction 
in the runs without surface fluxes in the clear region. Time-space plots of surface 
rainfall every three minutes are shown in Fig 2 and Fig 3 for TOGA COARE and 
GATE, respectively. The figures indicate that the precipitation patterns in the 
early stage of the storms are quite similar for all runs in each squall line case. 
However, the rainfall amounts started to diverge as the squall lines reached their 
mature stage (about 4 hours). The runs without any surface fluxes and without 
surface fluxes in the clear area showed a greater decrease in surface precipitation 
than the run with surface fluxes only in the clear area. 

Model simulated rainfall is also separated into stratiform and convective 
components according to the method proposed by Tao et al. (1993). Table 3 lists 
the stratiform rainfall amount and percentage of the total rainfall that is 
strtiform for all the runs. The results from TOGA COARE and GATE are not 
consistent. The TOGA COARE runs indicate that the surface fluxes increase 
stratiform rainfall, while the GATE runs do not show this. Instead, the 
percentage of the stratiform rainfall showed mixed results. This inconsistency 
was caused by small convective cells in front of the main squall line in the GATE 
runs (Fig. 3) in which the rainfall is mostly convective. Those secondary 
convective cells are caused by gravity wave from the main squall line convection 
which will be discussed in the next section. 

3.2 Cloud structure 

To examine the cloud structure and cloud coverage, hydrometeor contours 
are plotted for the TOGA COARE runs at 8 hours in Fig. 4 and the GATE runs at 
12 hours in Fig.5. Since the maximum updrafts are greater in the GATE case (12 
ms-1) than in the TOGA COARE case (8.5 ms-1) at this time (Fig. 6), the TOGA 
COARE squall line has much shallower convection (4.5 km) compared with the 
GATE squall line. The cloud coverage response to the surface fluxes is similar for 
both cases. Fig. 5 shows that cloud coverage is larger for the runs with surface 
fluxes (control run Gl) and with surface fluxes in the clear region (G3) compared 
to those runs without surface fluxes (G2) and with surface fluxes in the cloudy 
region only (G4). Cloud coverage in the TOGA COARE case (Fig. 4) is also slightly 
larger for the runs with surface fluxes (Tl) and with surface fluxes in the clear 
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region (T3). The increase in cloud coverage was caused by the moisture input 
from the sea surface evaporation in the cases with surface fluxes. Differences in 
cloud coverage for these various runs are consistent over the whole time period 
of simulation time. 

Although the squall line precipitation amounts responded similarly to the 
surface fluxes in similar trend, there was a difference in the convective cell 
generation ahead of the main line system between the TOGA COARE and GATE 
cases. The GATE case had stronger peak updrafts (Fig. 6) and taller convective 
cores compared with the TOGA COARE case (see Figs. 4, 5). The TOGA COARE 
simulations did not produce any convective cells ahead of the main squall 
system because of the weaker and more uniform updrafts, while the GATE 
simulations generated several convective cells in front of the main squall line 
(Fig. 3). These convective cells are generated by gravity waves excited from the 
convective cores of the squall line. As indicated in Fig. 6, the runs with surface 
fluxes (Gl) or with surface fluxes in the clear area (G3) have greater vertical 
velocities compared with runs G2 and G4. Therefore the Gl and G3 had more 
secondary convective cells ahead of the main squall line. 

3.3 Convective Available Potential Energy (CAPE) 

Fig. 7 and Fig. 8 show the space-time variation of the CAPE values in the 
TOGA COARE runs (Tl, T2, T3, T4) and the GATE runs (Gl, G2, G3, G4), 
respectively. The CAPE values were computed from model results at each half 
hour using GEMPAK (General Meteorological Package, desJardins et al., 1991). 
First, there are significant differences in the low CAPE (< 100 J/kg) region 
between the TOGA COARE and GATE runs. The GATE runs show much broader 
low CAPE bands and wake compared with the TOGA COARE runs. These 
differences are primarily caused by the more intense convection in the GATE 
case. Fig.9 shows the averaged profiles of temperature and moisture differences 
from the initial soundings averaged over the 150-km cloud band. The GATE case 
shows much greater mid-tropospheric warming and lower tropospheric drying, 
and therefore larger area of low CAPE after the convection. 

The main effects of the surface fluxes on the CAPE values are in the clear 
region. With environmental surface fluxes turned on in runs Tl and Gl and 
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runs T3 and G3, the CAPE values show less variation from start to the end 
compared to runs T2 and G2 and runs T4 and G4. Table 4 show that the CAPE 
values increase gradually with time from 1418 to 1586 J/kg for run Tl, while the 
CAPE decreases little from 1625 to 1350 for run Gl. This discrepancy between 
TOGA COARE and GATE is probably due to two reasons: the larger air-sea 
surface temperature difference and hence the larger surface fluxes in the TOGA 
COARE case and the stronger convergence of the GATE case. The runs without 
surface fluxes (T2, G2) in the clear region showed the largest decrease in CAPE in 
the clear region with respect to time due to the lack of moisture supply from the 
sea surface. 

Our analysis indicated that latent heat flux is the dominant source of CAPE 
compared to sensible heat flux. Four supplementary runs were carried out for 
this purpose. In the TOGA COARE case, the average value of CAPE in the clear 
region decreases from 1418 to 812 J/kg over the 16 hour simulation without 
latent heat flux, while CAPE increased slightly from 1416 to 1504 in a run 
without sensible heat flux. Similarly in the GATE case, a run without latent heat 
flux yields a large decrease in CAPE from 1625 to 850 J/kg, while a run without 
sensible heat flux had a relatively modest decrease in CAPE, from 1600 to 1300 
J/kg. This indicated that surface fluxes, especially latent heat fluxes are a major 
source of CAPE. 

It appears that CAPE values in the clear region are closely related to the total 
amount of surface precipitation in both the TOGA COARE and GATE cases. For 
example, in runs Gl and G3, the CAPE values are similar (Fig.7), and there is 
little difference (8%) in surface rainfall between those two runs. However, CAPE 
values decreased greatly from start to finish in runs G2 and G4, and the rainfall 
amounts had a greater percentage decrease (26%, see Table 3) as well. Similar 
conclusions apply to the TOGA COARE case. Since both the TOGA COARE and 
the GATE squall lines propagated quickly, the CAPE values in the clear regions 
out ahead of the systems determined the parcel buoyancy and vertical wind 
speed and therefore the intensity of the convection. 

3.4 Trajectory analysis and water budget analysis 
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As described in the previous sections, the precipitation, the domain 
averaged CAPE, and the cloud structure, all show remarkable similarity for the 
simulations with surface fluxes (control) and with surface fluxes in the clear 
region. This indicated that surface fluxes from the cloudy region had only a 
secondary influence on cloud development compared with those from the large 
clear area, even though the surface fluxes in the cloudy region can be almost 5 to 
7 times that of the clear area (see Fig 1). The cloudy area is too small to have a 
significant amount of moisture input. To support this conclusion, further 
analysis on the simulations has been performed. 

Fig. 10 and Fig. 11 show backward trajectories from the TOGA COARE and 
GATE control runs (T1 and Gl). The trajectory analysis used the model- 
simulated results from every three minutes and the data was linearly 
interpolated. The primary inflow of moist air is from in front of the squall line, 
and that air rises through the convective cores of the squall line. Nine parcels 
from the boundary layer out in front of the main squall line were randomly 
chosen. The trajectories indicated that cloud parcels mostly originated from the 
convective leading edge, and parcels from the leading convective edge are mostly 
from the clear area in front. These parcels ascended rapidly in a narrow 
convective updraft up to 4 to 6 km before getting transported to the stratiform 
region of the cloud system. 

Moisture enters the squall line convection mainly from two locations: 
advection of moist boundary layer air from out ahead of the squall line and 
surface moisture flux from directly beneath the convection. The boundary layer 
moisture in the clear area comes directly from evaporation of the ocean surface. 
Sensitivity tests have indicated that moisture flux is the most important factor in 
sustaining the CAPE. The amount of water vapor flowing into the convection 
from these two sources can be estimated from the cloud model results and from 
the surface flux calculations. Water vapor influx from the leading edge was 
estimated by integrating the horizontal moisture flux up to a 2 km height using 
the storm-relative U and the corresponding water vapor mixing ratio q in front 
of the squall line (about 4 km ahead of the convective core, see Fig.12). Table 5 
gives average moisture inflow values per unit meter of squall line computed 
from the three-minute cloud model data for both TOGA COARE and GATE. The 
moisture flux directly beneath the cloud contributed less than 8% of the total 



moisture. These results show why the surface flux in the large clear area has a 
greater effect on the squall line precipitation. More than 92% of the moisture 
comes from the boundary layer of the environment, which is ultimately from 
the surface moisture flux. The above analysis indicates that correct specification 
of fluxes in the clear region is very important in tropical squall line system 
modeling. 

4. Summary and Conclusions 

A two-dimensional cloud-resolving model is linked with a TOGA COARE 
flux algorithm to examine the impact of ocean surface fluxes on tropical squall 
line development and its associated precipitation processes. The distribution of 
CAPE in the model domain was also examined over the course of the cloud 
development. Although the cloud structure and precipitation intensity are 
different between the TOGA COARE and GATE squall line cases, the effects of 
the surface fluxes from different regions on the amount of rainfall and the cloud 
development processes are quite similar. The simulated total surface rainfall 
amount in the runs without surface fluxes is about 67% of the rainfall simulated 
with surface fluxes. The model results also indicate that the surface fluxes from 
the large clear environment are the dominant moisture source for tropical squall 
line development. Surface rainfall was reduced 8 to 9% in simulations without 
surface fluxes in the cloud region. Trajectory and water budget analysis indicated 
that most moisture (92%) is from the boundary layer of the clear environment 
air even though moisture fluxes are large in the cloudy region. Tropical squall 
lines feed on high CAPE air in the front of the convection, and they leave a low 
CAPE wake behind them. The conclusion is that accurate specification of the 
surface fluxes in the large environmental area is more critical when simulating 
air-sea interactions. 

In the tropical atmosphere, CAPE is relatively uniform and moderate over 
wide areas. Comparisons between soundings made in clear areas remote from 
cumulus clusters with those made in clear spaces within cloud clusters showed a 
steeper lapse rate in the ones far from cloudy areas (Bunker et al., 1949; Malkus, 
1958). Surface energy fluxes from the sea steepen the lapse rate, while transports 
by cumulus clouds work to restore a most adiabatic lapse rate (Riehl, 1954). 
Model simulations of tropical squall lines agree with those conclusions. As a 



continuation of this study, the 3D version of the GCEmodel will be used to 
simulate the PBL structure in both the clear and disturbed regions of various type 
of cloud systems (i.e. , fast/slow-moving squall system, and rotating MCS's) for 
comparison with CO ARE observations. Trajectory analysis of both the inflow 
and the outflow of various cloud systems along with changes in the ambient 
CAPE will also be documented. 
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Table Captions 


Table 1 Initial environmental conditions expressed in terms of CAPE, 
precipitable water, sea surface temperature (SST), surface air 
temperature, water vapor and wind for the TOGA COARE and GATE 
squall cases. 

Table 2 Setups for the eight sensitivity experiments conducted. 

Table 3 Surface rainfall amounts accumulated over 16 hours and normalized 
against a control run (T1 and Gl). The amount as well as the percentage 
of rainfall that was stratiform are also given. 

Table 4 The average CAPE (J/kg) in the clear area at the end of each sensitivity 
test. 

Table 5 Moisture contributions (kg s' 1 m' 1 ) from the cloudy region and from the 
leading edge of the squall system. 



Figure Captions 


Fig. 1 Spatial variation of (a) rain rate, (b) surface frictional velocity, (c) latent 
heat flux (LH), and (d) sensible heat flux (SH) at 10 hours into the 
simulation. 

Fig. 2 Space-time distribution of surface rain rate (mm/hr) for the TOGA 
COARE runs (Tl, T2, T3, T4). 

Fig. 3 Space-time distribution of surface rain rate (mm/hr) for the GATE runs 
(Gl, G2, G3, G4). 

Fig. 4 Vertical cross-sections of total hydrometeors at 8 hours into the TOGA 
COARE simulations (Tl, T2, T3, T4). 

Fig. 5 Vertical cross-sections of total hydrometeors at 12 hours into the GATE 
simulations (Gl, G2, G3, G4). 

Fig. 6 Maximum updrafts (every 3 minutes) in the TOGA COARE (Tl, T2, T3, 
T4) and GATE (Gl, G2, G3, G4). The traces were smoothed by a running 
average of 5 points. 

Fig. 7 Space-time plots of the CAPE (J/kg) in the TOGA COARE runs (Tl, T2, 
T3, T4). 

Fig. 8 Space-time plots of the CAPE (J/kg) in the GATE runs (Gl, G2, G3, G4). 

Fig. 9 Average temperature (DT) and moisture differences (DQ) from the 
initial soundings. The average was taken over the 150 km cloudy areas 
in runs Tl and Gl. 

Fig. 10 Backward trajectories from 9 parcels from the TOGA COARE control run 
(Tl). 

Fig. 11 Backward trajectories from 9 parcels from the GATE control run (Gl). 



Fig. 12 Schematic diagram showing the computation of moisture fluxes from 
front of the squall line and from the surface beneath the cloudy region. 



TABLE 1 



CAPE 

(J/kg) 

Precipitable 
water (g / cm 2 ) 

SST (°C), Tsfc (°C) 
Qsfc (g/kg), Usfc 
(m/s) 

TOGA 

COARE 

1418 

6.05 

28.0, 26.8 
20.01, 3.20 

GATE 

1625 

4.80 

26.9, 26.2 

17.33, 0.43 


TABLE 2 



TOGA COARE 

GATE 

Surface fluxes + 

LW + SW radiation (Control 
run) 

T1 

G1 

No Surface fluxes 

T2 

G2 

No surface fluxes in 
Cloudy Region 

T3 

G3 

No surface fluxes in 
Clear region 

T4 

G4 

















TABLE 3. 



TOGA CO ARE 

GA' 

IE 

■i 

Run number 

T1 

T2 

T3 

T4 

G1 

G2 

G3 

G4 

Domain average 
total rainfall 
(mm/1 6hr), 
Percentage 
relative to control 
run 

5.27 

(100%) 

3.53 

(67%) 

4.84 

(92%) 

3.78 

(74%) 

9.36 

(100%) 

6.30 

(67%) 

8.48 

(91%) 

6.97 

(74%) 

Stratiform rainfall 
(mm/ 1 6hr), 
Percentage relative 
to its total 

2.48 

(47%) 

1.27 

(36%) 

2.27 

(47%) 

1 

2.71 

(29%) 

2.33 

(37%) 

3.13 

(37%) 

1.81 

(26%) 


20 



























Table 4. 


TOGA COARE 
Initial CAPE: 1418 (J/kg) 

GATE 

Initial CAPE: 1625 (J/kg) 

T1 

T2 

T3 

T4 

G1 

G2 

G3 

G4 

1586 

754 

1510 

858 

1350 

810 

1302 

854 
















Table 5 



moisture from 

moisture from 


leading 

edge of squall line 

cloudy region 

GATE 

137.6 

9.8 

TOGACOARE 

70.4 

6.1 
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